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ABSTRACT 

We discover a definite correlation between the frequency of tlie quasi-periodic 
oscillations (QPO) in quiescence states and the duration of the quiescence state 
of the transient X-ray source GRS 1915+105. We find that while the QPO 
frequency can be explained with the oscillation of shocks in accretion flows, 
the switching of burst to quiescence states (and vice versa) and their duration 
can be explained by assuming an outflow from the post-shock region. The 
duration of the quiescence state is inversely related to the QPO- frequency. We 
derive this relation. We also flnd the correlation between the observed low 
(~ 0.001 — O.OlHz) and the intermediate (1 — lOHz) QPO frequencies. Our 
analytical solutions are verifled by analyzing several days of public-domain data 
from RXTE. 

Subject headings: accretion, accretion disks — black hole physics — 
hydrodynamics — shock waves — stars: individual (GRS 1915+105) — X rays: 
stars 
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1. Introduction 

X-ray transient source GRS 1915+105 in our galaxy exhibits various types of 
quasi-periodic oscillations with frequencies ranging from ~ 0.001 — O.OlHz to ~ 67Hz 
(Morgan et al, 1997; Paul et al. 1998; Yadav et al. 1999). The object is sometimes 
in a flaring state with regular and quasi-regular bursts and quiescences, while at some 
other time it is in usual low-hard and high-soft states. While the light curves look very 
chaotic with no apparent similarity between observations in two different days, some of 
the features are classiflable: (a) low-frequency QPO {u^ ~ 0.001 — O.OlHz) is due to the 
transition between burst and quiescence states (which we term as 'on'-state and 'off'-state 
respectively) and vice versa; (b) the intermediate frequency QPO (z// ~ 1 — lOHz) could 
be due to oscillations of shocks located at tens to hundreds of Schwarzschild radii Rg 
(= 2GM/c^ is the Schwarzschild radius. Here, M is the mass of the black hole, G, and c 
are the gravitational constant and velocity of light respectively) and (c) very high frequency 
QPO {uh ~ 67Hz), if at all present, could be due to oscillations of the shocks located at 
several Rg. vj is generally observed during quiescence states. Typically, a shock located at 
Rs (unless mentioned otherwise, measured hereafter in units of Rg)., produces an oscillation 
of frequency, 

1 1 D-aC^O ,-.. 

vi = — ~ —Rg — s (1) 

tff R Rg 

where, R is the compression ratio of the gas at the shock. Here we used a result of Molteni, 
Sponholz & Chakrabarti (1996, hereafter referred to as MSC96) which states that the 
time-period of QPO oscillation is comparable to the infall time (tm/a// = tffR ^ R^J"^) 
in the post-shock region. However, we assume now that the post-shock velocity is not 
necessarily R~^^'^ dependent as in a free fall but could be slowly varying, especially when 
angular momentum is high. In this case, tin fail oc i?" . Clearly, a ~ 3/2 for a low angular 
momentum freely falling matter and a ~ 1 for a post-shock flow of constant velocity Vqc/ R. 
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Here Vq is a dimensionless quantity which is exactly unity for a free-fall gas. For a gas of 
7 = 4/3, R ^ 7 and for 7 = 5/3, -R ~ 4, when the shock is strong. Thus, for instance, 
for a z// = 6Hz, Rs ~ 38 for M = IOMq and 7 = 4/3. For uh = 67Hz, i?, ~ 8 for the 
same parameters. MSC96 and Chakrabarti & Titarchuk (1995, hereafter CT95) postulated 
that since black hole QPOs show a large amount of photon flux variation, they cannot be 
explained simply by assuming some inhomogeinities, or perturbations in the flow. 

If the QPOs are really due to shock oscillations, they should almost disappear at 
low energy soft X-rays, since these X-rays are produced in pre-shock flow which does not 
participate in large-scale oscillations. Second, a shock-compressed gas with compression 
ratio R > 1, must produce outflows or extended corona which pass through sonic points 
located at at Re = foRs/2, where /o = R^/{R — 1), if the flow is assumed to be isothermal 
till Re (Chakrabarti 1998, 1999, hereafter C98 and C99 respectively). In this solution the 
location of the sonic point Re and ratio between outflow and inflow rates are functions of 
the compression ratio R of the shock alone. Till the sonic point Re, matter is subsonic and 
this subsonic volume is fllled in a time of (C99), 

AttRI < p> 



Hill 



- - (2) 



where, < p > is the average density of the sonic sphere, and M^y^t is the outflow rate. The 
Compton cooling becomes catastrophic when < p > Rekes ~ 1, ^es = 0.4 is the Thomson 
scattering opacity. Thus the duration of the off-state (i.e., duration between the end of a 
burst and the beginning of the next burst) is given by, 

t.U = ^=^. (3) 

We use now a simple relation between inflow and outflow rates given by (C98, C99), 

Mout n ©o-ut-Rr R i3/2 /3 R . , . 
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where, 6j„ and Qout are the solid angles of the inflow and the outflow respectively. Because 
of the uncertainties in 9j„, Qout and M(„ (subscript 'in' refers to the accretion rate) we 
deflne a dimensionless parameter, 

yJm iVlEdd 

where, MEdd is the Eddington rate. Using Eqs. (4-5), we get the following expression for 

toff as, 

_ 10.47 RlR]exp{h-l) 

Or, eliminating shock location Rs using eq. (1) and a = 3/2, fo = 1, we obtain. 



^off 



expifo-i) , M 



V3„-4/3, 



For an average shock of strength 2.5 ^ R ^ 3.3, the result is insensitive to the compression 
ratio. Using average value oi R = 2.9 and for Qj^^ ~ 0.1 (which corresponds to 0.1 
Eddington rate for Qout ~ ©m) we get. 

Thus, the duration of the off-state must go down rapidly as the QPO frequency increases if 
the flow geometry and the net accretion rate remains flxed. When one considers a constant 
velocity post-shock flow, a = 1, and vq = 0.066 (chosen so as to keep the same numerical 
coefficient as in eq. 8) the above equation is changed to, 

t = 461.5(^)(^^)-^(^^)2z/72.. (9) 

^^ ^Qm IOMq' ^0.066^ ^ ^ ' 

Interestingly, Vq = 0.066 (i.e., a constant velocity of seven percent of the velocity of light) is 
very reasonable for a black hole accretion. 

If the hot post-shock gas of height ~ Rg intercepts n soft photons per second, from the 
pre-shock Keplerian component (CT95), it should intercept about nf^/A soft photons per 



-6- 

second when the sonic sphere of size Re is filled in. Thus, the photon flux in the burst state 
should be about /q /4 -^ 4 times larger compared to the photon flux in off-state. Depending 
on the degree of flaring, and the fact that the wind is bent backward due to centrifugal 
force, the interception may be higher. 

Since ul is basically due to recurrences of on and off-states, it is clear that 
z/i ~ 1/ < toff + ton >■ Here, ton is the duration of the burst state which may be very small 
for extremely regular (spiky) bursts reported by in Taam, Chen & Swank (1997) and Yadav 
et al. 1999. In this case, 

., = 0.0022(|^)(i^).,^ Hz. (10) 

When on-state has a non-negligible duration (ton 7^ 0), it is found to be directly related to 
the toff (Belloni et al., 1997; Yadav et al. 1999). Assuming ton ~ ^o//, the z/^ would be less 
by a factor of two when on states are broad. When the burst is very regular but 'spiky' 
(i.e., with momentary on-state), ton ~ 0. The presence of ul for these regular 'spiky' bursts 
are reported in Manickam and Chakrabarti (1999a) 

Thus, if our shock oscillation solution for QPO is correct, the observations must pass 
all the following tests: (a) the QPO in the off-state must disappear at low energies, (b) the 
QPO must generally be absent in the on-state, when the sonic sphere is cooled down, (c) 
the intermediate QPO frequency must be correlated with to// as in Eqs. (8-9) and (d) the 
photon flux must jump at least a factor of 4 or more when going from quiescence to burst 
state. In addition, (e) lowest frequency u^ observed must be correlated to the intermediate 
QPO frequency uj by eq. (10). There are uncertainties regarding the inflow velocity and 
actual volume-flUing time, but we expect that above relations to be satisfied in general. 

In the present Letter we show that observations do pass through these tests and 
therefore the shock oscillation model may be the correct picture. In the next Section, we 
present detailed analysis of some of the observational results on GRS 19 15-1-105 available in 
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public archive and show how they point to the shock oscillation model. Finally in §3, we 
make concluding remarks. 

2. Observational Results 

Figure 1 shows a light curve of the first phase of observation of June 18th, 1997, on 
the right panel. The average count rate (per second) vary from around ~ 5000 in the 
off-state to about ~ 24, 000 in the on-state. The ratio of the fluxes is about ~ 5. The 
duration of these states vary chaotically. At the mean location of a few off-states (arrows 
on right axis), observation time is marked in seconds. For each of these off-states, the 
power density spectrum (PDS) in arbitrary units is drawn in the left panel. The most 
prominent QPO frequencies {uj in our notation) are connected by a dashed curve just to 
indicate its variation with time. There are some weaker peaks which follow the short dashed 
curve, indicating that they may be higher harmonics. Observations of this kind for several 
other days show similar variations in QPO frequencies and details are presented elsewhere 
(Manickam & Chakrabarti, 1999ab). 

In Figure 2, variation of uj with the duration toff of the off-states (triangles) for the 
whole observation period on June 18th, 1997 in the log-log scale. Observational results 
from several other days (May 26th, 1997; June 9th, 1997; June 25, 1997; October 7th, 
1996; October 25th, 1996) are also plotted on the same curve with circles, filled squares, 
squares, filled circles and stars respectively. We did not put error bars since in duration 
scale is it uniformly ±2 seconds, and in frequency scale error bar is decided by the chosen 
bin-size while obtaining the PDS (In our analysis it remains around 0.15-0.25 Hz, increasing 
monotonically with QPO frequency). Times at which the photon flux is halved during the 
on-to-off and off-to-on transitions are taken respectively to be the beginning and the end of 
an off-state. Equation (9) is plotted in dashed lines with (from uppermost to the lowermost 



line) Qm = 0.0034, 0.0123, 0.0163, 0.028, 0.0293 and 0.043 respectively indicating a slow 
variation of the accretion rate, provided the collimation property remains the same. Two 
dotted curves, on the other hand, represent eq. (8), with 6^ = 0.06 (top) and 0.093 
(bottom) respectively. We find that the inverse-squared law (eq. 9) may be a better fit to 
the observations. Since on Oct. 7th, 1996 the points are lumped closed to the lower right 
corner, not much could be said about whether it follows our relation or not, but it is to be 
noted that its general behaviour (low frequency, high duration) follows our result for any 
reasonable Q^- 

Table 1 shows the variation of the uj (taken from Fig. 2) with days of observations. In 
3rd column, the expected ul has been put {vl/'^ for Oct. 7 and June 18 results as they 
show ton ~ ^o//)- In 4th Column, the observed Vl is given. Generally, what we observe is 
that, when the drift of z// is large, PDS around vl is also broad In any case, observed vl 
agrees with our expectations. 

Figure 3 shows the PDS of the off-state centered at 1576s (see Fig. 1) of the June 18th 
observation. The energy range is given in each panel Clearly, QPO disappears completely 
at low energies, exactly as is expected in the shock oscillation model (MSC96) though QPO 
frequencies, when present, seem to be energy independent (see, Manickam & Chakrabarti, 
1999b for details). The pre-shock fiow which emits soft radiation participates little in 
the oscillation as the fractional change in the Keplerian disk due to shock oscillation is 
negligible. On the contrary, the fractional change in the size of the post-shock flow during 
the oscillation is large. Thus, the flux of hard X-rays oscillates as the size of the post-shock 
region oscillates. This is the cause of QPO in our model. In on-states, when they exist, the 
QPO is found to be very weak. 
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3. Discussion and Conclusions 

In this Letter, we have discovered a relation between the QPO frequency in 1-lOHz 
range with the duration of the quiescence state at which the QPO is observed. We also 
derived a relation between the low QPO frequency and the intermediate QPO frequency. 
We analyzed several days of RXTE observations and showed that our relations are satisfied, 
especially when the average bulk velocity in the post-shock region is constant. We showed 
that the QPO disappears in the low energy, but is very strong in high energies. The photon 
flux is found to fluctuate, typically by a factor of 4 or more, indicating that a vertically 
inflated post-shock region is responsible for interception of the soft photons from a Keplerian 
disk. This factor seems to be similar to /q /4 (C98, C99) for any reasonable compression of 
the gas, which strengthens our belief that the quasi-periodic cooling of the sonic sphere of 
the outflow from the post-shock region may be responsible for the rapid transitions between 
on and off-states. Our computation of the duration of the off-states from this considerations 
are found to be quite reasonable. We find that for z/^^/^ law, toff is insensitive to the mass 
of the black hole while for z/~^ law, toff is inversely proportional to the mass. Trudolyubov 
et al. (1999) recently found the duration of 'hard' states varies as —7/3 power of the lowest 
centroid frequency for a group of data while we find an inverse-squared law when we choose 
the QPO frequency where the power is strongest. 

Although we chose a specific model for the outflow (C99) for concreteness, the physical 
processes invoked are generic and the explanation should be valid even when other models 
for outflows are used (except self-similar models). The shock location near the inner edge 
of the Keplerian disk can drift on viscous time scale (see appendix of CT95 where the 
transition from Keplerian to sub-Keplerian is plotted as a function of viscosity). The shocks 
can evacuate the disk, and form once again, very similar to what was seen in the numerical 
simulation of Ryu et al. (1997). This drift would cause a drift in frequency as Trudolyubov 
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et al. (1999) recently showed (see also, Belloni et al. 1997; Markwardt, Swank & Taam 
1999, Muno et al. 1999). Our model invokes also outflows (which we believe form naturally 
in the post-shock region and from the transition region from Keplerian to a sub-Keplerian 
flow) which we find useful to explain variations in photon counts between 'off' and 'on' 
states. 

This work is partly supported by a project (Quasi Periodic Oscillations in Black Hole 
Candidates) funded by Indian Space Research Organization (ISRO). The authors thank 
NASA for making RXTE data available and ISRO for creating a Data Bank at their Centre 
where these data are stored. 
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TABLE 1 
Correlation Between Low And Intermediate Frequencies of QPO (in Hz) 

Day Vi Ul from Eq. 10 Vl observed Remarks 

Oct. 7, 1996 1.7-1.90 0.0013-0.0016 

Oct. 25, 1996 2.9-3.5 0.0006-0.0009 

May 26, 1997 5.32-6.54 0.0076-0.0115 

June 9 1997 6.22-7.42 0.0138-0.0197 

June 18 1997 4.0-7.51 0.003-0.011 

June 25 1997 6.61-7.42 0.016-0.02 



0.001 


broad on 




sparse data 


0.01 


spiky 


0.014 


spiky 


0.01 


broad on 


0.022 


spiky 



- 12- 

REFERENCES 

Belloni, T, Mendez, King, A.R., Van der Klis, M. & Van Paradijs, J. 1997, APJ 488, L109 

Chakrabarti, S.K. & Titarchuk, L.G. 1995, ApJ 455, 623 (CT95) 

Chakrabarti, S.K. 1998 in Observational Evidence for Black Holes in the Universe, ed. S K 
Chakrabarti 19 (Dordrecht: Kluwer Academic Publishers) (C98) 

Chakrabarti, S.K., 1999, A&A, 351, 185, (C99) 

Manickam, S.G. & Chakrabarti, S.K. 1999a in Proceedings of Young Astrophysicists of 
Today's India, Ind. J. Phys. 73(6), 967. 

Manickam, S.G. & Chakrabarti, S.K. 1999b A&A (submitted) 

Markwardt, C.B., Swank, J.H. & Taam, R.E. 1999, ApJ 513, L37 

Molteni, D., Sponholz, H. & Chakrabarti, S.K. 1996, ApJ 457, 805 (MSC96) 

Morgan, E. H., Remillard, R. A. & Greiner, J. 1997, ApJ 482, 993 

Muno, M.P., Morgan, E.H. & Remillard, R.A. 1999, ApJ (submitted) 

Paul, B., et al. 1998, ApJ 492, L63 

Ryu, D., Chakrabarti & Molteni, D. 1997, ApJ 474, 378 

Taam, R.E., Chen, X.-M., Swank, J.H., ApJ, 485, L83 

Trudolyubov, E. Churazov & M. Gilfanov, A&A Letters 351, L15 

Yadav, J.S., et al. 1999, ApJ 517, 935 



This manuscript was prepared with the AAS L^T[t;X macros v4.0. 



-13- 



Fig. 1. — Plot of the light curve (right panel) and evolution of power density spectrum (left 
panel) of the first phase of June 18th, 1997 observation. Off-states analyzed are marked 
by the time of observations on the right axis. The QPO frequencies (where the power is 
strongest) are connected by a dashed curve to highlight the evolution of z// with time. 



Fig. 2. — Variation of QPO frequency vj with duration of quiescence states to//- Dotted 
curves are the to// oc v~*/^ law (eq. 9) derived using simple free-fall velocity assumption. 
Dashed curves are the toff oc z/^^ (eq. 10) with constant post-shock velocity law. The general 
agreement strongly points to the shock oscillation model. 



Fig. 3. — Power density spectrum of the off-state (duration 87s) centered at 1576s constructed 
from selected channel intervals of the binned PGA data. QPO is seen only in high energies, 
strongly pointing to the shock oscillation model. 
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